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Based on a set of cosmological N-body simulations we analyze properties of the dark matter haloes (DM) 
in a galaxy mass range (10 11 — 10 13 h -1 M@) in modified ACDM cosmology with additional dynamically 
screened scalar interactions in DM sector. Our simulations show that scalar interactions support picture 
of the Island Universe. Rapid structure formation processes are shifted into higher redshifts resulting in a 
much smaller accretion and merging rates for galactic haloes at low redshifts. Finally, we present how this 
"fifth" force affects halo properties, like density profile, triaxiality, ellipticities and the spin parameter. 
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1 Introduction 

In this paper we study the impact of the long-range scalar interactions (LRSI) in the dark matter (DM) 
sector 11121 on some properties of the dark matter haloes. LRSI was proposed to overcome some difficulties 
that ACDM model faces at scales important to the galaxy formation. Large-scale clustering of LRSI 
cosmologies was studied in (3J-0. Here we will study DM haloes density profiles, mass accretion histories, 
shape and spin parameters. 

2 Theory & numerical application 

We model the LRSI in DM sector by introducing modifications to the Newtonian gravity. Effective gravi- 
tational potential mimicking fifth force between two DM particles of mass m is 

$(r) = -^(l + /3e-^) , (1) 

Here G is Newton's constant; r and x = r/a(t) are, respectively, the particle separation vector in real 
and comoving coordinates. The j3 parameter measures the relative strength of the scalar-interaction com- 
pared to the usual gravitational interaction, while the r s parameters is Yukawa-like screening length and it 
measures the effective range of the scalar interactions. Accordingly, the modified potential gives rise to a 
modified force law between DM particles. The new expression for the force law is 
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To study the structure formation and collapse of DM halos we use publicly available GADGET2 code 
by Volker Springel J6). GADGET2 is a hybrid N-body code that uses two different methods of calculating 
the particle forces. At large scales it uses a Particle-Mesh approach, an algorithm in which the Fast Fourier 



* Corresponding author E-mail: pchela@camk.edu.pl. Phone: +48 223296 102, Fax: +48228410046 

Copyright line will be provided by the publisher 



2 



W. A. Hellwing: Galactic halos in cosm. with LRSI 



Transform method is used on a uniform lattice grid to solve the Poisson equation for gravitational potential. 
Then the potential is used to derive particles accelerations. We use Eqn. (fl]i to modify the Green's functions 
in the k space to obtain modified gravitational potentail. At small scales the code use Oct-Tree algorithm for 
obtaining the particle forces. We have implemented the Eqn. (f5]) in GADGET2 to obtain the modified force- 
law (see Ref. J4] for the details). Finally to find haloes and calculate their properties from our simulations 
snapshots we have used the Amiga Halo Finder (AHF), kindly provided to the public by A. Knebe & 
S. Knollmann J7J. 



3 Results 

We have preformed a series of N-body simulations in boxes of width 25h~ 1 Mpc containing 256 3 DM 
particles. Therefore the mass resolution (mass of a single particle) in our numerical experiments was 
m p ~ 7.7 x 10 7 /i _1 M Q . Cosmology used in our simulation assume fi m = 0.3, ft a = 0.7, Clk = 0, h = 
0.71 and cr 8 = 0.8. 

To allow for a direct comparison between corresponding halo in different models, we identify a galactic 
halo in the ACDM run and then we cross-corelate (c-c) this halo with another halo found in the simulations 
with LRSI. As a criteria for matching halo we use: (a) - the same origin in space (in vicinity of 5 x r V i r at 
redshift of first identification); (b) - proximity of the final space occupation; (c) - adjacency of the velocity 
vector at z = 0; (d) - the mass criteria (related halo in the LRSI model cannot excess five times the original 
ACDM halo mass). 

We begin by studying the shape of the density profiles of the c-c haloes. In the Fig. Q] on the left 
we show density profiles (solid lines) of the c-c haloes together with the best fits to the NFW profile [8] 
(dotted line): Pnfw = Ps/[{r/r s )(l + r/r s ) 2 ], where p s and r s are scaling density and scaling radius 
respectively. Arrows mark the distance from the center of the halo corresponding to the doubled force 
resolution 2 • e ~ 20/i _1 kpc. We decided to not trust the results, below this distance, because they 
can be highly affected by two-body scattering. We note that for that particular halo NFW fit is not so 
good, especially for a halo in the model with f3 = 0.2, r s = l/i _1 Mpc, and (3 = 0.5, r s = l/i -1 Mpc. 
Nevertheless it is striking to depict, that the modified force-law and the modified potential do not change 
the universal shape of the halo density profile, known from the ACDM N-body simulations. 
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Fig. 1 On the left: Density profiles of the cross-corelated haloes (solid lines). Dotted lines indicate NFW best fit. 
Arrows show the distance corresponding to the two times the force resolution « 20h~ 1 kpc. On the right: Mass 
accretion histories of the cross-corelated haloes (solid lines). Dotted lines coresponds to the best fit of the equation 
<[3j) to the data. Horizontal dashed line indicates the final mass at z — 0. In each panel at the bottom-right corner we 
present the formation redshift Zf derived from best fit to the data. 
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Table 1 Summary of the properties of the cross-corelated halo in a different models at z = 0. M is the mass of halo 
\10 12 hT 1 Mq\, a v is the 3D velocity dispersion [kms^ 1 ], v max is the maximum of the rotation curve \kms~ 1 ], r V i r 
stands for the virial radii [h^ 1 kpc], T is the triaxialites parameters, ei and are ellipticities, A is the spin parameter 
and Zf is the formation redshift as defined by the Eqn. {3]l. 



Model 








M 


o v 


^max 




T 


ei 


e-2 


A 


*/ 


ACDM 








4.18 


258 


247 


328 


0.65 


0.24 


0.15 


0.044 


0.9 


P = 0.2, 


r s 


= lh- 1 


Mpc 


5.12 


303 


266 


351 


0.93 


0.23 


0.21 


0.047 


1.95 


P = 0.5, 


r a 


= lh- 1 


Mpc 


6.96 


389 


322 


389 


0.76 


0.28 


0.2 


0.048 


2.66 


P = 1.0, 


r a 


= 0.5/1" 


_1 Mpc 


6.86 


454 


312 


387 


0.41 


0.37 


0.13 


0.103 


1.18 


p = 1.0, 


r s 


= lh- 1 


Mpc 


8.34 


517 


388 


413 


0.51 


0.12 


0.06 


0.032 


3.67 


p = 1.0, 


r s 


= 2h~ l 


Mpc 


12.15 


572 


426 


468 


0.3 


0.13 


0.037 


0.053 


3.82 



On the right side of the Fig. Q] we plot mass accretion histories (MAH's) of the c-c haloes (solid 
lines) versus the expansion factor a. Following the procedure presented in (9) we fit MAH to the simple 
exponential law, given by: 

i i 2 

M(a) = exp [a(l - a )] , z f = 1 , (3) 

where M = Mq is the halo mass at z = 0, a is the scaling parameter, a is the cosmic scale factor, 
and Zf is the formation redshift defined with the help of the a parameter. We also plot best fits to MAH's 
for a single c-c halo. From these plots we see that c-c haloes in cosmologies with LRSI tend to have much 
higher formation redshifts Zf. Except cases with p = 1, r s = 500/i~ 1 kpc and /3 = 0.2, r s = lft^Mpc, 
c-c haloes in modified cosmologies have much quieter recent accretion histories then the ACDM case. 
Especially when we look at MAHs of c-c haloes with P = 1 and ?' s = 1 or 2/i _1 Mpc, we note that these 
haloes have not experienced any major merger since a ~ 0.4 (z = 1.5), and then they accreted mass 
steadily until present time. 

Some properties of the cross-corelated haloes are given in the Table Q] Note that corresponding halo 
in presence of the scalar interactions gets more massive at the final redshift. The value of the velocity 
dispersion (a v ) or the virial radiu£] r v are also higher compared to the ACDM case. Now we define halo's 
shape and spin parameters as: 

a 2 - b 2 c b J 
triaxiality: T=— -; ellipticities: e\ — 1 e 2 = 1 ; and spin: A = 



a a y ZAl v i r v v i r r v i r 

where we adopted definitions of the shape parameters from IflOl and the spin parameter from IfTTI . Here 
a > b > c are eigenvalues of the inertia tensor of a halo, J is the total angular momentum of the halo 
and M V i r ,v V i r and r V i r are its virial mass, velocity and radius respectively. Inspection of the properties 
collected in the Tab. [T]indicate that there is a big scatter between c-c haloes of their shape parameters, but 
we can note a weak tendency of higher A for stronger scalar interactions (except the case for P = l,r s = 
l/i^Mpc). 

Now we turn to the statistical properties of the haloes formed in our simulations. We average three 
shape parameters and the spin parameter over all haloes found in simulation boxes which satisfy this mass 
criterion M > lO 11 /! -1 A/©. We present results in Table|2] Clearly, the presence of the scalar interactions 
does not change shapes of the haloes. Triaxiality and ellipticities are converging to the values of the ACDM 
case within one sigma errors. On the other hand we can note, that halos in cosmology with the LRSI, tend 
to have higher spin parameters. 



We define virial radius as the radius at which the density of the halo drops to the value A ~ 340 times the background density. 
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Table 2 Mean ellipticities, triaxiality and the spin parameter when averaging over haloes with M > 10 11 /! Mq. 



Model 








< ei > 


< e 2 > 


< T > 


< A > 


ACDM 








0.25 ±0.09 


0.14 ± 0.09 


0.58 ± 0.22 


0.040 ±0.024 


P = 0.2, 


r s 


= Ih- 1 


Mpc 


0.25 ±0.08 


0.15 ±0.08 


0.59 ±0.22 


0.044 ±0.028 


13 = 0.5, 


r s 


= Ih- 1 


Mpc 


0.25 ±0.08 


0.14 ±0.09 


0.57 ±0.23 


0.047 ±0.029 


13 = 1.0, 


r s 


= 0.5/i " 


_1 Mpc 


0.26 ±0.08 


0.16 ±0.09 


0.60 ±0.21 


0.052 ±0.03 


13 = 1.0, 


r s 


= lh~ 1 


Mpc 


0.25 ±0.09 


0.14 ±0.08 


0.57 ±0.21 


0.054 ±0.032 


13 = 1.0, 


r s 


= 2/j" 1 


Mpc 


0.24 ±0.08 


0.13 ±0.08 


0.55 ±0.23 


0.053 ±0.04 



4 Summary 

We have performed N-body simulations of structure formation in cosmologies with the additional long- 
range scalar interactions in the dark matter. We have studied the impact of such modification of the standard 
ACDM cosmological model on the c-c haloes and on the averaged statistical properties of the halo family. 
We conclude that the outer parts of the c-c halo density profile do not show any significant deviations from 
the ACDM halo density profile. We have to note though, that the c-c halos of our choice had their virial 
radii smaller then the screening length of scalar interactions models considered here. It will be interesting 
to study the inner most parts of the density profiles of cross-corelated haloes. To do that a much higher 
resolutions in forces and in masses are needed. We have also studied the MAH's of the c-c haloes. Here we 
see clear imprint of the LRSI on the mass accretion histories. Haloes in modified gravity have much higher 
formation redshifts, and they tend to show less mergers at small redshifts. This indicates that, in cosmology 
with the LRSI, small-scale structures like DM haloes form at earlier epochs and support the picture of 
galaxies being island universes at low redshifts. These results are encouraging in the view of apparent 
difficulties of reconciling the predictions from the ACDM cosmology galaxy formation simulations with 
astronomical observations (see Ref. JSJHIEI f° r detailed discussion). We have also measured averaged 
shape and spin parameters of the halo populations. We conclude that LRSI does not affect the shape 
of the haloes. However we see a clear signal in the spin parameter. DM haloes with LRSI have more 
angular momenta than their cousins in the ACDM case. The research presented here is just a first step in 
understanding how additional long-range interactions in the dark matter can affect various properties of the 
gravitationally bounded DM haloes. To study this matter in greater detail a more advanced study is needed. 
We have showed that modified gravity imprints on average halo spin and on particular halo mass accretion 
histories. Therefore it will be, in principle, possible to distinguish between the ACDM and LRSI models 
using data from future detailed galaxy surveys. 
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